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Abstract 

It is shown that gauge field-dependent fermion Dirac operators 
from lattice QCD form an ergodic operator family in the probabilistic 
sense, provided the gauge field is an ergodic random field. As a con- 
sequence, the integrated density of states of such Dirac operators in 
the thermodynamic limit exists and is almost surely independent of 
the chosen gauge field configuration. 
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1 Introduction 

The Standard Model of Elementary Particles is one of the biggest achieve- 
ments of theoretical physics during the second half of the 20 th century. It 
provides a common conceptual basis for all elementary forces except gravity. 
The part which describes the strong nuclear force is called Quantum Chromo- 
dynamics (QCD). The associated Lagrangian density has a clear and simple 
appearance that fits on a single line, 

C(x) = ~F ia ,{x)F' w {x) + ^(x)(iD(A)-m)i/j(x), (1) 

with D(A) := 7^(9^ + iAy) being the Dirac operator of the fermion field 
i/j(x), 4>(x), which depends on the gauge field A ■= (A^x))^, and F^ u (x) = 
d^A u (x) - d u A^(x) - [A tl (x),A u (x)] being the field tensor. In spite of its 
structural simplicity concrete quantitative predictions are difficult to derive 
from ((T]), and one often resorts either to calculations in perturbation theory 
or numerical simulations on a discretized (Euclidean, after Wick rotation) 
space-time, known as lattice QCD (LQCD, see [H] for an overview). 

Several basic properties of QCD, such as the spontaneous breaking of 
chiral symmetry or the phenomenon of quark confinement, manifest them- 
selves in the regime of low energies, where perturbation theory in the QCD 
coupling constant cannot be applied. Spontaneous chiral symmetry breaking 
is signaled by the formation of a non- vanishing chiral condensate In 
a seminal paper [TJ, Banks and Casher formulated a link between the value 
of the condensate and the spectral properties of the Dirac operator D in 
the deep infra-red. Since the gauge field A does not appear explicitly in the 
observable, it acts as a background field which nonetheless determines the 
spectral properties in a non-trivial way. 

The idea that the distribution of the low-lying eigenvalues of the fermion 
Dirac operator is very close to the one of the corresponding (i.e., respecting 
symmetries) random matrix ensemble was put forward in [TTJ [161 E] an d 
affirmed by numerous numerical studies, e.g. [61 [2], for a review see for 
example [TJ)]. In fact, these distributions agree to an accuracy that would, 
perhaps, allow to replace the derivation of average spectral properties of 
the fermion Dirac operator by sampling gauge field configurations with the 
random matrix eigenvalue distribution. The robustness of this phenomenon 
over a broad range of parameter values, like the underlying gauge group 
or the system's temperature (in the Boltzmann weight) is also remarkable 
and encourages us to formulate our ultimate goal as to establish a rigorous 
mathematical link between these distributions beyond numerical evidence. 

The present paper is a first modest step towards this goal, namely to 
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fit the formulation of the model in LQCD into the mathematical framework 
of ergodic operator families which has been built up over the past three 
decades or so to study random Schrodinger operators and, especially, the 
Anderson model. In contrast to random Schrodinger operators, however, the 
randomness in LQCD models lies on the lattice bonds - not on the lattice 
sites - and corresponds to a random magnetic field, rather than an alloy or 
a quenched glass. 

As the main result of this paper we prove that Dirac operators of LQCD 
which depend on the gauge field indeed constitute ergodic operator families 
in the probabilistic sense, provided the gauge field itself is ergodic, i.e., has 
sufficient rapidly decaying correlations (see Section 11.1.21 for a precise for- 
mulation). This, in turn, is a fair assumption in many physical situations, 
e.g., at high temperature or if the gauge field is massive (which is believed 
to be true for non-abelian gauge fields). As a consequence of our result the 
integrated density of states exists in the thermodynamic limit and is almost 
surely independent of the chosen gauge field configuration. 

Many observables can be expressed in terms of derivatives of the QCD 
partition function with respect to source terms. For example, the chiral 
condensate is given by [15J 

(W) = -lim lim ^-d m \ogZ QCD , (2) 

m->0 V->°o y 

with the partition function 

Z QCD = r Y\ dA^ det[iD(A)+m]e- SYM( ^. (3) 

J xeV,n=l,...,d 

Here, the integration of the fermionic variables yields the fermion determi- 
nant det[iD(A) +m], and Sym(A) is the Euclidean Yang-Mills action. 

It is customary to use the quenched approximation in numerical simula- 
tions, which amounts to setting the fermion determinant is equal to one. This 
reduces the numerical effort significantly and corresponds to the physical case 
of infinitely heavy sea quarks. 

The discretization of the Dirac operator is also subtle, because the naive 
discretization leads to the occurrence of fermion doublers, which have no 
physical meaning. There are several ways to work around this problem. 
Wilson proposed to add a term that vanishes in the continuum limit and 
suppresses the doublers on the lattice [T7j. Another method is to introduce 
staggered fermions - the lattice is divided up in sub-lattices where different 
staggered phases live, that are interpreted as physical phases [H H3]. We 
are mainly interested in those two cases, where the Dirac operator still has 
nearest-neighbour interaction. This is not the case for another elegant solu- 
tion, the overlap operator proposed in [9]. 
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1.1 Introduction of the Model 

Now we come to the precise description of the mathematical setting. We 
consider the lattice Z d , d>2, with fermion fields supported on the sites and 
gauge fields supported on the bonds of the lattice. We take only the one 
particle case into account. Then the matter fields are complex vectors and 
the configuration of all matter fields is supposed to be an element of the 
Hilbert space - £ 2 (Z d ,C k ) of square summable Z d -sequences in C k . 
is equipped with the usual scalar product 

M) = E I>(z)^(z), <p,ipeJf- (4) 

X<=Z d 1=1 

We will consider operators on Jtf, that also depend on the configuration 
of gauge fields on the bonds of Z d . The set of bonds in the lattice Z d is 
denoted by 

B ■= Z d x{l,...,d}. (5) 

The bond (x, //) e B is the one connecting x and x + e M , with e M the unit 
vector in Z d pointing in direction fi. We give the bond (x, fj,) the orientation 
from x to x + e M . 

The gauge fields associated to the bonds are elements of a compact Lie 
group G, the gauge group. We assume Q to be either SO(N), SU(N), or 
U(N) to be explicit and since these are the relevant physical cases. The 
gauge field on the bond b = (x,/i) is denoted by U x>fl or XJ h . 

It turns out to be necessary to consider the direction of a bond, the gauge 
field for going from x + e M to x is U x+ ^,-n and we set U^,-^ = U X) ^ 1 . 

A gauge field configuration is the collection {Ub}^B- As mentioned before, 
this gauge field configuration is randomly generated. To specify the under- 
lying probability space we will need the notion of a plaquette, a collection of 
four bonds that form a plane square in Z d , 

P(x;n,u) := {(x,ii),(x + e^v),(x + e 1/ ,iJ,),(x,v)}, (6) 

with fj, + v. We need the product of the gauge fields along a plaquette 
p = p(x;n,v), 

Up ■= U X ;ix,v '■- U x u U x+ £ v ^ U x+ e^ t u U Xjfl , (7) 

where the orientation of the bonds leads to the inverse gauge fields. Thus we 
have U x -^v = U~} v and define a plaquette as positively orientated if \x < v. 
The set of all positively orientated plaquettes is denoted by 



V ■= {p(x; 12,1s) | x e Z d , /j,, v e {1, . . . , d], [i<v}. 



(8) 
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1.1.1 The Probability Space 

We start by specifying the probability space for a single bond. Since Q is a 
compact Lie group, Q is equipped with a natural measure, the Haar measure 
(jlh- Because of the compactness of Q, the Haar measure [in is normalized 
and we obtain the probability space (G,J-i,(J<h), with T\ being the cr-algebra 
of Borel sets of the topological space Q. A priori, we let the gauge field U b on 
a single bond be a random variable that is uniformly distributed with respect 
to the Haar measure fin of Q. 

The gauge field configuration U = {U b } b ^B can be regarded as an element 
of the product space Q B . Q B is compact in the product topology by Ty- 
chonov's Theorem. The probability space for gauge field configurations is 
now constructed by means of cylinder sets as in [3J. A cylinder set is subset 
of Q B of the form 

M = {Ueg B \U bl eA 1 ,...,U bn eA n } (9) 

with Ai, . . . , A n 6 T\ and b\, . . . , b n e B. The set of all cylinder sets is denoted 
Z. We take as a a-algebra for Q B the cr-algebra T generated by the system 
of all cylinder sets. The probability-measure P on Q B is then defined to be 
the product measure, setting for every cylinder set 

n 

¥({Ueg B \U bl eA l ,...,U bn eA n }) = YiMA)- (10) 

8=1 

Now, we modify the measure P by a weight function that represents the 
gauge action. Formally, the measure P is defined as 

dF(U) = Z- X e~ s{ ^ df(U) (11) 

with Z being a normalization factor. We assume the gauge action to be of 
the following form: 

S(U) = /3£ReTr(l-Lg (12) 

with f3 > and U p the plaquette variable as defined in ([7]). This action on 
discrete space-time is the Wilson action used in lattice QCD calculations. 
Note that the gauge action is invariant under translations. Denoting by T e 
the translation in Z d by £ e Z d , i.e. T t x := x - £, and defining the translation 
of a gauge field by I e Z d to be 

T l U x ^ ■■= U T i x ^ = U x -i ilM , (13) 



for any U X)fX e Q and (x, fi) e B, we have 

S(U) = S(T e U) 



(14) 
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for all gauge field configurations U_^Q B and all I e Z d . 

In the following, we use the Gibbs formalism to establish the existence of P 
and its uniqueness for small (3, following |3]. To this end, we fix a finite subset 
A of B and a gauge field configuration rj e Q B that represents the boundary 
condition, which fixes the gauge field outside of A. The corresponding local 
specifications are 



d^XpiH) = {Zip) X exp 
where 



-j8 E ReTr(l-[/ p ) 

pnA±0 



Y\dMUb) (15) 

feeA 



Y[dMU b ). (16) 

6eA 



/3 X! ReTr(l-£/" p ) 

pnA*0 

It is known (see for example [121 H]) that, for sufficiently small /3 > 0, these 
measures weakly converge, as A / B, to a unique Gibbs measure, since C/ B is 
compact, the action is continuous in the gauge fields, and for all b e B there 
is a constant d < oo, such that 

£ HReTrtl-tgL < d. (17) 

Furthermore, the Gibbs measure is independent of the boundary condition 
rj. In our case, we have 

£ HReTrCl-tgL < 2(d-l)-2J\T, (18) 

peV-bep 

uniformly in b € B, since any bond 6 e 25 is element of 2(d - 1) plaquettes. 

Dobrushin's uniqueness criterion ensures the uniqueness of the Gibbs 
measure for translation invariant interactions, provided 

< £ (|p|-l)||ReTr(l-£/ p )L < /T 1 , (19) 

peV:bo€p 

for one bond (and, hence, for all bonds) b e see for example [T2]. In our 
case, the criterion is fulfilled for all 

< P < n r- (20) 

Thus the Gibbs measure P exists and is unique, for sufficiently small j3 > 0. 

We note in passing, that the translations are measure preserving trans- 
formations with respect to P and P, that means for all A e J 7 , I € Z d 

¥(T i A)=F(A) and P(TM) = P(A). (21) 

Put differently, P and P are stationary w.r.t. the group Z d of translations. 
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1.1.2 Ergodic Probability Measures 

A stationary probability measure P is called ergodic iff, for all A, A' e J 7 , 

—1— £ P(AnT l A')^P(A)P(A'), asL-oo, (22) 

^ + ^ ' ZeZ", \\l\\^L 

with || Z|[ co = max{|/i|, . . . , \l d \}. 

A random variable / : (G B ,F) -> (M.,B) is called invariant iff f(T e U_) = 
f(U_), for all £ e lt d and almost all U_ e Q B . The importance of the notion 
of ergodicity lies in the fact that any invariant random variable is P-almost 
surely constant. 

In the following, we show that the measure P is ergodic, provided (120]) 
holds true. First, we specify the decay of correlations. We define a metric 
d ■ B x B -» RJ on B by setting d(b, b) := and 

d(6,6') := (23) 
min {n | 3{p 1} . . . ,p n } c P : 6epi, pi n p 2 * 0, • • • ,Pn-i n Pn * 0, 

for b,b' £ B, b ± U . I.e. <i(6, 6') is the minimal number of plaquettes to connect 
b and b' . 

\{p = (x; fi, u), p = (x; jl, u), and pnp + then fx-x 1^ < 1. Therefore, the 
minimal number of plaquettes connecting b = (x,/jl) e B and b' = (y,^) e B is 
at least fx -y\\oo- Observing that, for all x e Z d , /i, i/, r e {1, . . . , d}, v + n + r, 

p(x; /i, z/) n p(x + e^; fi, v) + 
p(x; fx, v) np(x - e M ;/x, v) + 

p(x;/jL, v) r\p(x;fj,,r) + 0, (24) 

we obtain, with b = (x, /i), b' = (y, u) e B as above, that 

||a;-y||oo < d(b,b') < \\x-y\\i + d, (25) 

with \\x -y\\i = Eti - 

A modified version of the metric d called d ■= ln(-|)<i is obtained by multi- 
plying d with ln(f ), where c> (3. 

Now, we take two cylinder sets A, A' e T and show that condition (122]) 
is fulfilled for A and A'. Since the sigma- algebra T is generated by Z, (T22]) 
extends to all T. There are two finite sets Aa,A A / c B and A b ,A' b , e T\ for 
all b € A a, fr 7 6 A A' such that 



(26) 
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Now choose XA,e, XA',e e C(Q B , [0, 1]) such that XA,e> XA',6 depend only on the 
variables Iff, with b e A A or 6 e A A ', respectively, XA,e(L0 = 1? XA',e(U') = 1; 
for all £^,£'6 A', and 

F{U e A c | X A,eOD > 0} < e , P{C/' e | X ^, £ (£') > 0} < e. (27) 

Then XA,e,XA',e are continuous functions that differ from the characteristic 
function of A, A' only on an set of measure less than e. 

Now, we use a result of [5] summarized in [12] . It states that, if 7 < 1, 
where 7 is a constant depending on the interaction, one gets for any two 
bonds i,j e B that 

f XA,, XA',e dF - f X A,e d¥ f X A>,e dl 
Jg B Jg B Jg B 

< \ e-^)(l -7)- 1 \(xA,e) *j(xa'J, 



where 



and 



A,(/):=£e^)sup{ 1/00-/(^01 I ^ = C/ 6, & (29) 



ieB 



7 = sup X! e<^'W (30) 
In our case for z j, can be estimated as 

Pij< £ ||/5ReTr(l-[/ p )||<2iV/3 1[^,j) = l], (31) 

such that we get 7 < 3 • 2(d - l)2A^/3 • | ,and 7 < 1 corresponds to 

c < l2iV^T)' (32) 

arriving at Condition f )20|) . We denote the distance of ^4 and ^4' by 

dist(A,A') :=min{(i(z, e A A ,j e A A ,} (33) 
and the diameter of A by 

D A :=max{d(i,j)\i,j e A A }. (34) 

If z € A A then 

A,( XAe ) < j;^)i[ fc6Ail ] < |A A |(|)^ (35) 

keB P 
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and analogously Aj(xA',e) ^ l-^-A'Kf )'° A ') provided j e A A i. Inserting this into 
Equation (128|) . and choosing i e and j e A A > such that d(i,j) = dist(A, A'), 
we estimate 

f XA,e XA', e d¥ - f X A,e dF f X A>,e dF 
Jg B Jg B Jg B 

, x dist(A,A')-(D A +D A ,) I a I I a | 

- llcj i-uN(d-iy (36) 

for all < (3 < c< (12N(d - l)) -1 . In the limit e 0, we obtain 

(„x dist(X,A') 

c) ' (37) 

where 

C ^'4 1-127V( d -l) c ' A ^'b) < 38 ' 

is independent of the distance of A and A'. 

The exponential decay of correlations implies at once, 

1 Y, \F(AnT e A')-¥(A)¥(A')\ 

etz d , \\e\\oo<L 

i „x dist(A,T £ A') 



(2L + l) d 



< 



(2L + iy £eZdj i^i^sl 



m-dist( J 4,A')-2(DA+-D^') 



, 2^//3V dlSt(A ' A ' ) " (DA+ ^ ) 1 

- ~2ltt[7) t^^ ' (39) 

where we used that d(k,T e k) <d(i,k) + d(i,j) + d(T £ k,j) and thus 

dist(A,T l A') = mm{d(i,j)\teA A ,j eT £ A A ,} 

> mm{d(k,T e k)\k e A A r} - max{<i(z, k)\i e A A , k e A^} 
-max{d(T £ A;, j)\k e A A >,j e T^A A '} 

> ||£|U - - 2Z?^ - dist(A, A')- (40) 

The validity of (I2"2"j) on cylinder sets extends to the sigma-algebra T gen- 
erated by the cylinder sets Z by a monotone class argument. Thus Condition 
(1201 for the uniqueness of P ensures its ergodicity, too. 
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1.1.3 Ergodic Families of Wilson Dirac Operators 

We specify the considered operators. The dependence of those operators Djj 
on the gauge field configuration is emphasized by the index U_. We consider 
the corresponding family of operators {Du_} Ue gB. 

Let {Du} U£ gB be a family of bounded, self-adjoint operators on Jif = 
£ 2 (Z d ;C k ). We call this family stationary if it depends on the gauge field 
configuration U_ € Q B in such a way that translations act transitively, i.e., 

T e DuT- £ = D T tu (41) 

for all U_^Q B and £ € Z d , where r l denotes the corresponding translation on 
J^ 7 , i.e., 

[t £ </>](x) = <Kx-1), (42) 

for any <fi e x e Z d . A stationary family {Djj_} Ue gB is called ergodic if 
the underlying probability measure P on Q B is stationary and ergodic. The 
crucial fact about ergodic families {Du_} Ue gB is the independence of their 
spectra on U_, P-almost surely. [10] . 

We assume that Dij_ includes only nearest-neighbour interaction, i.e., for 
cf) € J^, x € Z d , [Du<f)](x) depends only on the values of <p(x), <f)(y) for those 
y with \x - y\ = 1 and the gauge fields U x ^, U x -^ for fj, € {1, ... , d}, where we 
use the notation V x -n '■= £4-e M ,At- 

There are various examples for such operators of physical interest. As 
mentioned in the introduction, we are mainly interested in the Wilson Dirac 
operator and the staggered fermions operator. For simplicity, we concentrate 
our attention to the Wilson Dirac operator in this paper. 

In lattice gauge theories the Wilson Dirac operator D is used [8], which 
is a discretized version of the QCD-Dirac operator Dqcd - 1^(9^ + iA^) + m 
with gauge fields A^. The corresponding matter fields are defined on the 
hypercubic lattice Z 4 and are assumed to have a Dirac structure labeled 
by Dirac indices a e {1,2,3,4}, as well as a colour structure with labels 
c e {1, . . . , N c }. The Dirac structure is represented by the 4x4 Euclidean 
Dirac matrices {7^)^=1,. ..,4. A customary explicit representation is 





7.= * ^ (43) 



-201,2,3 

I 

with o"i.2,3 being the Pauli matrices. The Dirac matrices form a Clifford- 
Algebra since they fulfill {7^,7,,} = 25^. Introducing 75 := 71727374, i-e., 
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we observe that {7^,75} = 0. 

The gauge group is Q = SU(N C ) and acting on the colour structure. 
Therefore, k = 4 • iV c , and (/) = {0 a , c }a=i r ..,4, c =i...,jv c e £ 2 (Z 4 ,C k )- The Wilson 
Dirac operator D is defined by 

[D0] QjC (x) := Y, A x ) 
0=1 { 

4 N c 1 

" K E E E (( r 75)a,/3-^(7 5 7^)a,/3) (^ W )c,/ 00,/(2: + O^) f, (45) 

/i=l(7=±l/=l J 

in short, 

4 

[ZtyKx) = 7 5 [0(a:) - « £ E ( r " ^T/J^w <t>(% + o^)] . (46) 

The parameter r e (0, 1] is the Wilson parameter and k > the hopping 
parameter. 

Displaying the dependence of D on the gauge field configuration U_ by 
writing Djj, we observe that Djj fulfills condition ( l4"Tj) for any {7 e SXJ(N C ) B 
and any € Jjf, x e Z d , 

[r^r-V](^) = [%r-V](x-0 

4 

= [t - V] E( r + ^)^-^ [t' 1 4>]{x- i + ae^) 

[1=1 cr=±l 

4 

= (j)(x) -KY, E ( r + 1<rn)U x -t,<Tii 4>( x + ™p) 

fl=l <T=±1 

= [iW>](x) (47) 
and hence we have 

r^ZW = £W, (48) 

for all £ € Z d and £ e S[/(iV c ) B . Thus, if P is ergodic, SO is {-D[/}[7e5(7(Ar c ) s ) 
and its spectrum is P-almost surely constant. 

1.2 Introduction of the Integrated Density of States 

In the following we study the integrated density of states of {-Dj/jt/sgB, which 
represents the number of eigenstates per unit volume. For the precise defi- 
nition of the integrated density of states, we restrict our analysis to a finite 
subset A c Z d . Besides, this also allows us to relate our analysis to numerical 
simulation. 
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The boundary of A, denoted <9A, is defined as 

«9A := {yeA\3xeA c , \x-y\ = l} c A. (49) 

Furthermore we use the canonical orthonormal basis ■= {s^ x '^} x€ z d ,i€{i,...,fc} 
of Jt? where the Z d -sequence s^'*) is set to be 

^•' , (//): ('" (50) 
[0, x + y 

with ej e C k the unit vector in direction i. 

Since D(y contains only nearest-neighbour hopping, the value of Du_s^ x ^ 
does not change, if we replace Djj by a restriction of Djj_ to A, for any point 
x in A\<9A. Only the boundary dA needs further specification. We present 
two customary choices for this, namely, Dirichlet boundary conditions and 
periodic boundary conditions. 



1.2.1 Dirichlet Boundary Conditions 

First, we restrict the operators to the finite subset A c Z d by means of the 
projection 

Pf r) :Jf^ A W , [P[ dir) v](x):={ X£A ' (51) 

10, x i A, 

with 

^ dir) = £ 2 (A,C k )cJf (52) 
being the Hilbert space of sequences vanishing outside A. Note that 

P TlA = r e P A r- e . (53) 

Then we define 

D ua = P A ir) DuP[ dir) ■■ J? A (dir) - M> K {dir) • (54) 

Note that, since Jif A dir ^ is finite-dimensional, D^ l A can be represented 
by a matrix of size (&|A|) x (/c|A|), where |A| denotes the number of elements 
in A. Since Djj_ is self-adjoint, so is D^f 1 ^ . The number of eigenvalues of 
Dy 11 ^ smaller than some E e R, counting multiplicity, is denoted by 

N A f\E) := {![£><*> <E]}. (55) 
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The integrated density of states of D^f^ is defined as the number of eigen- 
values smaller than E per unit volume, 

Pf2^) := ±N^(E). (56) 

Clearly, N A dl j\E) depends only on the gauge fields on the bonds connecting 
points in A. Note that probabilistic statements about N^J^ U (E) do not 
depend on £ eZ d , since P and Djj_ are stationary. 

1.2.2 Periodic Boundary Conditions 

Another way to restrict Djj_ to a finite set A c Z d is to require periodic 
boundary conditions, which is often used in numerical simulations. In order 
to define periodic boundary conditions we assume A to be a cube of side 
length L. Without loss of generality we may assume that A = {1, 2, ... , L} d . 

We define the Wilson Dirac operator on £ 2 (A( per ); C fc ) with periodic boun- 
dary conditions by 

[D ( £ r U](x) := 7 5 [>(*)-^E E (r-^)^w (Kx + aeJ] , (57) 

where := (Z/LZ) d and x + ae^ is determined only modulo multiples of 

L in all directions. Similarly 

= (58) 

where x - e M is also defined modulo L. Thus D^ er ^ : J^ per - ) -> J^ per \ with 

J? A (per) := f(A^;C k ), (59) 

only depends on the values of Ufor bonds 6eAx{l,...,d}, i.e., on 

U A ■= {U x ,„} X eA,^i,..., d . (60) 

Since Jtf A per ^ is finite-dimensional, we can transcribe the definition of the in- 
tegrated density of states to periodic boundary conditions. We set N A P ^\E) 

to be the number of eigenvalues, counting multiplicity, of D^ e A smaller than 
E e ffi. and define the integrated density of states in the periodic case as 

P ( £L\E)--=^N^\E). (61) 
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2 Main Theorem 

Our aim is the definition of the integrated density of states for {Du_}u_ e ^. A 
natural way is to let A be a cube of side length L and investigate the case 
L -> oo, the thermodynamic limit. As it turns out, the boundary conditions 
imposed is immaterial. 

Theorem 2.1. Let (Q 13 ,^, P) be the probability space defined in Section [7T71 
and choose < (3 < 12 jv(d-i) ( su °h that P is ergodic). Let {Du} U£ gB be the 
family of Wilson Dime operators on Jtf. Let (Q n ) n m be a sequence in Z d of 
nested cubes, fl n £ fl n +i, with Q n / Z d . 



(i) Then the limits 



and 



^(£0 == iim^<:2(^) (63) 



n-»oo 



exist for all i?eR, F-almost surely, and are independent of the sequence 

(ii) Furthermore, for all the integrated density of states p(E), de- 

fined by 

p^\E) = P ^\E) =: p(E), (64) 

zs independent of the chosen boundary condition and of U_ e P- 
almost surely. 



3 Proof 

3.1 An Estimate on Eigenvalues 

Suppose, we take two disjoint sets f2i,f2 2 c Z d and an operator Djj_, that 
fulfills the requirements of Theorem 12.11 We can restrict Djj_ to f2 2 and 
Q\ u fl 2 as in (154)) by means of the projections P^, Pq 2 and Pn lU n 2 - Then 
we can determine the number of eigenvalues below some £el for all three 
restrictions, denoted by Njj^^E), Njj^^E) and Nu,q iU q 2 (E), respectively. 
If dist(Qi,Q 2 ) > 2, we know that Njj^ iU q 2 (E) = N EjUl (E) + Nu,u 2 (E), since 
the operator Djj_ links only neighbouring sites. Our first goal is the derivation 
of an upper bound on the difference of Njj^ (E)+Nu,q 2 (E) and Nu,a 1 un 2 (E). 
To this end, we start with a general observation for finite matrices. 
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Let A, B be complex, self-adjoint (M x M)-matrices. The rank of B 
is denoted by b, and the interesting case is b « M. Since A and B are 
self-adjoint, so is A + B, and all three matrices A, B and A + B have M 
real eigenvalues counting multiplicity. Due to the fact that rank(i3) = b, B 
has (M - b) eigenvalues equal to zero, and b eigenvalues different from zero. 
Furthermore, we denote by Na e No the number of negative eigenvalues of A, 
by Nb, N^ b , Na+b the number of negative eigenvalues of B, -B, and A + B, 
respectively. 

Lemma 3.1. Let A, B be self-adjoint M x M -matrices. The difference of the 
number Na of negative eigenvalues of A and the number Na+b of negative 
eigenvalues of A + B is at most rank(5) ; 

\N A -N A+B \<Tank(B). (65) 

Note, that the bound (165]) is independent of ||S||. 



Proof. First, we show that Na+b - Na < rank(S). Let us assume that 
Na+b > Na + rank(S). Then the min-max principle ensures the existence 
of a subspace X £ C A/ , with dimension dim(X) = Na + rank(S) + 1 such that 



sup {(j>\{A + B)(j>) <0. (66) 

cj>zX,\\<t>\\=l 



In particular we have 



sup {cf)\(A + B)<p) = sup {<p\A(p} <0. (67) 

</>eXnker(S), ||0||=1 4>eXnkeT(B),\\<j>\\=l 

Using the min-max principle again, we obtain Na > dim(X n kei(B)) > 
dim(X) - rank(5) = Na + 1. Therefore we have that Na+b ~ Na < rank(S). 

Now, we set A' := A + B, B' := B and get analogously N A '+b' - N' A < 
rank (5') that is Na - N a +b < rank (B). □ 

Lemma 3.2. Let (G B ,T, P) be the probability space defined in section Yl.l.l\ 

and choose < f3 < 12 jv(rf-i) suc h that P is ergodic. Let {Djj_} Ue gB be the 
family of Wilson Dirac operators on ffl . Furthermore let Qi, . . . ,Qj c Z d be 
disjoint, finite sets and Q := U/=i their union. 



(i) Then we have, for any U_e Q B and all E e 



1 

Pf 



<- * ^ m 
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(ii) If in addition the sets Qi, . .. ,Qj are cubes, such that Q is also a cube, 
then 



1 

Pf 



|0| 



(69) 



3=1 ~ 

for any U_eQ B and all E e R. 

Proof. We remark that it is enough to prove the case E = 0, because we 
can replace Djj_ by Djj_ - E. We start with the case of Dirichlet boundary 
conditions and define (k\Q\ x A;|f2|)-matrices A and C by 



(dir) 



(dir) 



(70) 



The matrices are chosen in such a way that we get with counting multiplicity 

<?(0)=Tr{l[A>0]} (71) 

and 



(72) 



E<S(0)=Tr{l[C>0]}. 

3=1 ~ 

Now, we set the matrix B := A - C to be the difference of A and C. The rank 
of 5 can be estimated as follows 



rank(£) < &£|<9%|, 



using that S = -fhj^-fb r By Lemma IBTTj we obtain 



1 

W\ 



3=1 



< k 



|n| 



(73) 



(74) 



and (i) is proven. 

To prove (ii), we use that, for any cube A, we have that 

rank [D™ - D$] < k \dA\ 

Therefore, (i) and another application of Lemma [3.11 yield ( 169 



(75) 



1 

Pi 



i / 



<- ¥ \{\n { £Z\e)-n^\e)\ 



N u,n ( E ) " E N uZj ( E ) 

3=1 3=1 



< 3k 



(76) 
□ 
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Lemma EOl is an estimate on the change of the integrated density of states 
as the subset of Z d is broken up into smaller pieces. The estimate is, indeed, 
precise enough to prove the existence of a limit in the sense of Theorem 12.11 
as is done in the next sections. 

3.2 Existence of the Integrated Density of States for a 
Special Sequence 

In this section it is shown that a limit for the integrated density of states 
exists almost surely for a sequence of growing cubes in Z d . 

To this end, we define the following sequence of growing cubes, 

A n :={-/ 2 n - 1 + l,...,/ 2 n - 1 } d , (77) 

with l e N to be fixed later. Note that A n has side length l 2 n . 

The virtue of the sequence (A n ) ne N is, that A n+ i splits into 2 d disjoint 
cubes, each of size |A n |, in a natural way. More precisely, there are z%, . . . , z 2 d e 
Z d such that, for 

n„,:={T^,...,T^} (78) 
being the set of associated translations, 

A n+1 = (J TA n . (79) 

TeU n 

In order to clarify the notation, we also introduce the sets U l n for / > n 
that consist of the translations needed to compose A; of translations of A n , 

: = {T n . . . T ; _i : T n € n n , . . . , e LVi}. (80) 

Thus, U n = II^ +1 and we have 

Al = (J TA n, (81) 

see Figure [0 Next, we study the integrated density of states of A n , as n 
grows. We omit the dependence of N^ r ^(E) and N^ r jj{E) on E and the 
gauge field configuration U_ and write 

N^[A n ] := N^l(E) and N^[A n ] := N^^E) (82) 

instead. 
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Am 



A 









Figure 1: A m and A n in Z 2 , with m = n + 2. The dotted lines indicate the 
translates of A n whose union gives A m . 

Lemma 3.3. For any l e N, the sequences (|^A^( dir )[A n ]) neN and 
(pQ A r ( per )[A„])^ eN converge, F-almost surely. 

lim t^M*0 [An] - ^ , lim t^JV^ [A„] - p<f r) . (83) 
Furthermore 

lim — I N^ dir) [AJ - N (per) [AJ 1 = 0, (84) 

n^oo |A n | 

IP- almost surely, and pi ■= = pj^^ ■ 

Proof. We show that (rhN^r) [A n ]) nefq is a Cauchy sequence and in the 
proof we denote N( dir )[A n ] =: iV[A n ]. The proof for periodic boundary con- 
ditions is completely analogous. 

Assume that m > n. By applying (IHTj) . one can split A m into 2 d ^ m -^ 
cubes of size |A„|, 

A m = [J TA n . (85) 



Ten: 



The mean integrated density of states for these translations of A n is 

E uV [TA " ] = £ ^ [TA " ] - (86) 

" TeII m r^-nl | J1 -m| Tell" 1 
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A 


An 

















Figure 2: Both sets A n and A m are split up in smaller cubes of the same size 
as A no . For A m only part of the splitting is sketched. 



Thus we can estimate 



1 N[A m ] - ^iV[A„] 



|A r 



|A r 



N[A m ]- Y, N[TA n ] 



J7) 



Ten" 



|A 



Y N[TA n ] - rr-MK] 



ni \ Tell" 



|A„| 



Lemma 13.21 directly gives us an upper bound for the first term on the right 
side of (1571). since we have 



|A r 



N[A m ]- Y N[TA n ] 



, h nd(m-n) 2(j(/ 2 n ) d 1 



< k 
2dk 



2 d ( m - n )|aA n | 



o 



(88) 



independently of the gauge field configuration. 

The second term on the right side of (1871) is the difference of the integrated 
density of states for a cube A n and its spatial mean over 

2 d(m-n) translated 

disjoint cubes of the same size. As we do not know, yet, whether this term is 
small with high probability, provided n is large enough, we split A n and its 
translates into smaller cubes of size |A no | for some no < n € N, as indicated in 
Figure [2J 
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We estimate 



|A 



' Y N[TA n ] - ^iV[A„] 



< 



m\ Tell" 1 

n 

1 



Ar 



1 



N[TA n ]- Y, N[T'TA no ] 



T'ell" 



ml Ten" 



2 iV[TA no ]-— J) N[TA no ] 



n\ Ten",, 



1 



Y N[TA no ]-N[A n ] 



Tell" 



39) 



using that any T e II™ is given as a product T = T'T", for unique T 7 € IT™ 
and T" e Il™ . Then Lemma I3~2l yields again an upper bound for the first and 
the third term on the right side of (|89|) . and analogously to (188|) . we obtain 
that 

' ^ N[TA n ] - Y N[T'TA no ] 

T'en" 



|A 



I 



m| Ten" 



2dk „ 
< 2~ n ° 



(90) 



and 



|A„| 



Y N[TA no ]-N[A n ] 

Tent? 



2c/A; „ 
< 2~ n °. 



(91) 



Thus Equations rtBB]). fjB5]) . (M)]). and 1151]) yield 



|A 



1 iV[A m ] - r^-MK] 



I A. I 



(92) 



'o 



1 



|A 



Y N[TA no ] 



m| Ten" 



|A 



"I Ten" 



We can choose no and then n > no so large that —(2 n + 2 n °) is arbitrarily 
small. To estimate the remaining term, we view {Z x (U)} xiZ d, with 



Z X (U) := N[T^A no ] = N^l oxAno ^E), 



(93) 



to be an invariant family of random variables. By Birkhoff's Ergodic Theo- 
rem, the mean of these random variables converges P-almost surely. Hence, 



m| Ten 



Y N[TA no ] 



(94) 



m=no+l 
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is a Cauchy sequence, P- almost surely. 

As noted above, we can replace iV[#] by N^ per \jf\ and repeat the proof 
for periodic boundary conditions with exactly the same arguments, since all 
sets are cubes and Lemma l3~27 ii) applies. 

Equation (|84|) is similarly proven as Lemma 13.2( h). 

J-\N^[A n ]-N^[A n ]\ < - 0, n-oo. (95) 

□ 

Note that, while the preceding lemma holds for all Iq £ No, this does not 
imply the independence of the integrated density of states of the choice of Iq. 
It turns out, however, that not only the independence holds true, but that 
furthermore the size of the cubes in the sequence is immaterial, as long as it 
is monotonically growing. 



3.3 Proof of Main Theorem 123 



The proof is similar to the one of Lemma 13.31 We choose Iq , n e N arbitrary, 
but fixed. Given Q k , there is an m e N such that Q k £ A m . We define 



S fc := {Ten™ |TA no cfi fc }, 
Q k := (J TA no c Q k 



(96) 



Note that Cl k is a rectangular box, whose smallest side length is at most 
two times smaller than its largest side length and all side lengths are multiples 
of /o2 n °. Moreover 



\n k \-\n k \ = \Q k \n k \<\dQ k \-\A no \ 

Now, we estimate for k > n , Q k 2 A no 



(97) 



|A 



' N[A k ]--^N[Q k ] 



|A 



' ^[ A *]-TTT E N i TA «o\ 



|A 



fc l Ten* 



|A 



1 



i 



£ iV[TA no ] - -~-N[n k ] 



-~-N[n k ] - r^-Mn,] 



(9* 
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In the first and third term we can apply Lemma 13.21 directly to get upper 
bounds. The second term converges P-almost surely to zero, by Birkhoff 's er- 
godic theorem. For the last term we estimate the difference of the integrated 
density of states of Q k and Q k , 



\\n k \-N[n k ] 


-p k 


\-N[n k ]\ 


p k 


■ p k \ 





p k \ \u k \ 
First, we estimate 

p k \-N[ti k ]-\n k \-N[n k ] 

\n k \(N[ti k ] - N[n k ]) + (p k \ - p k \)N[n k ] 

< \Q k \\N[Q k ] - N[Q k ]\ + (p k \ - |O fc |)iV[^]. 
Then we observe that 

\N[n k ]-N[n k ]\ < k\dn k \ + (2d+i)k(\n k \-\n k \) 



holds true because of the following estimate using Lemma 13.21 
\N[Q k ]-N[Q k ]\- £ N ^ 



N[Q k ]-iN[Q k ]+ N[{x}] 



< k(\dQ k \ + 2d(\Q k \-\Q k \)). 



Thus we get 



p k \ 



< 1 |iv[O fc ]-iV[Q jfc ]| + J_iV[fi fc ] 



p k 



p k \ 



(p k 


-p k \) 







< k\on k \ | (2d + 2)k(p k \-\n k \) n^oo Q 



(99) 



(100) 
(101) 



(102) 



(103) 



by using Equations fllOOl) . fllOip and the fact that p^iV[f2 fc ] < k. Altogether 
we have proven 



lim 



K\ N[Ak] -\k\ mk] 



= o, 



(104) 
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P-almost surely. The periodic case is again proven completely analogously. 
Recall that Lemma 13.31 gives the existence of the limit 

PuJE) ■- Km jijiOB) = «5l ]k\ N '*'Z (E) (105) 
with (A n ) neN as in (1771) . P-almost surely. Equation (11041) ensures the indepen- 



dence of this limit of the chosen sequence (£lk)ktN- Furthermore, Lemma 13731 
implies that the limit for periodic boundary conditions is the same. 

Since the integrated density of states is invariant under translations, it is 
P-almost surely constant, and (ii) follows. This finishes the proof of Theorem 

[ZD □ 
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